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Abstract The work is devoted to the implementation of the hydrodynamic laws to the head-on heavy ion collisions 


within the energy range 50-100 MeV/A. The hydrodynamic mechanisms of the bubble and ring structures formation 


are investigated. It is shown that there is a possible hydrodynamic explanation of the different structures being formed 


in the case of soft (K=200 MeV) and stiff (K=400 MeV) equations of state. Within the suggested approach the final 


geometry of the system is defined in the initial stage of the collision and is very dependent on the sound velocity in the 


nuclear matter. The obtained results are in a good correspondence with the Boltzmann-like transport theory 


calculations and the experimental data for the selected energy range. 
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1 Introduction 


The knowledge of the nuclear equation of state (EOS) 
is one of the fundamental goals in nuclear physics! ^l 
which has not yet been achieved. The possibility of 
extracting information on the EOS at high baryon 
density, is restricted to two fields of research that are 
observations of astrophysical compact objects and 
studies of hot nuclear systems, created in high energy 
proton induced reactions or in heavy-ion collisions 
(HIC). Probably this explains the intensive studies of 
heavy ion collisions during the last several decades. 
The very exciting parts of them are the studies devoted 
to the nuclear multifragmentation (MF) when a 
number of intermediate mass fragments 3<Z<20 (IMF) 
are being formed from the compressed and excited 


791 IMFs may provide a unique probe to study 


system 
the reaction mechanism and hot nuclear matter 
properties". In spite of such an interest to the topic, 


there are still many questions to be answered. The 


serious challenge in the studies is the inclusiveness of 
the heavy ion collisions experiments! !!. The nowadays 
detectors are very sophisticated and allow the high 
detection rate but still some data are missing! l. At this 
step a powerful instrument of computer simulations is 
of high importance. There is a number of transport 


[3-5]. and 


theories used for simulating collisions 
statistical decay models! that provide interesting 
microscopic results and give predictions on possible 
MF mechanisms. 

Within the Boltzmann-Uehling-Uhlenbeck 
(BUU) model it was shown that at the energy range 
60—75 MeV/A “bubble” or “doughnut” structures can 
be formed depending on the stiffness of the nuclear 
EOS'^3 [ater the 


structures was 


formation of “doughnut” 


experimentally! "?!. 


confirmed 
Unfortunately up to now there is no model describing 
in details the underlying physical mechanisms of the 
exotic shapes formation in the head-on HIC. There are 
different possible ways to look for the explanation. 


One is to continue with the computer simulations. In 
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that case one should bear in mind that we are still far 
from having models that are formally well founded, 
practically applicable, and sufficiently realistic to be 
quantitatively useful*!, Another way is to consider the 
and nucleon-nucleon 
look for the 
lor hydrodynamics 


qualitative Van der Waals 


Dll and to 


[22-23 


interactions similarity 
explanation in thermodynamics 
well developed for ordinary liquids. In that case a 
number of possible physical mechanisms being 
responsible for the phenomena can be suggested. 
Among them are the mechanism that includes the 
shockwave formation and interference and the one that 
studies the thermodynamic properties of the system 
(temperature, surface tension) together with the 
influence of the angular momentum on the system 
behavior". The other possibilities are the possible 
capillary effects with *capons""?! being formed on the 
nuclei surface and the spinodal decomposition of the 
system. Therefore, some deep analysis of the 
phenomena is needed to choose the appropriate decay 
channel. 

The first step toward a better understanding of 
the process is the analysis of the qualitative picture 
that comes from the BUU calculations and the 


[7-19] First to mention is the fact 


available experiments 
that the expanding velocity of the outer surface of the 
system is much smaller than the expanding velocity of 
the inner surface that suggests the importance of the 
surface effects. Second is the system behavior during 
the final BUU calculations 


simultaneous breakup with few fragments of similar 


stage. predict the 
masses and low kinetic energies. This suggests that the 
dynamic process might be responsible for the system 
evolution but at the breakup stage Coulomb force 
seems to be of great importance. Third is the angular 
distribution of fragments. In the case of soft EOS it is 
almost isotropic when for the stiff EOS fragments are 
in the plane perpendicular to the beam direction. 

From all the 
approach seems to be quite attractive as it allows 


above the hydrodynamic 


explanation of the observed behavior of the system. In 
this work, we report an attempt to develop the 


hydrodynamic shockwave 


[17,26] 


approach with the 
mechanism that can explain the observed 


phenomena and reveal its physical nature. 


2 Model description 


When studying the system from the hydrodynamic 
point of view the question arises whether such a 
macroscopic theory works for the nuclear systems. 
Originally the idea to apply such a macroscopic theory 
1970s-1980s by 
The domain of their 


introduced in 
[28] 


to nuclei was 
[27] 
and Stocker 


studies was in a higher energy range than we are 


Siemens 


interested in, nevertheless, the existing confirmation 
for the applicability of hydrodynamics for nuclear 
systems suggests that it should be possible to use it for 
the energy range in focus"?! 

The symmetry of the system allows the 
simplification of the model. Namely from the collision 
of two nuclei it is possible to change to the collision of 
a spherical nucleus of radius R with a rigid wall that 
moves toward it with the velocity vo. The difference is 
in the viscous effects, therefore, the slip boundary 
condition on the wall surface should be added. 

In our studies we start from the acoustic 
approximation developed in the theory of liquid 


031] Within this model according to 


droplet collisions 
Huygens principle the expanding nucleus edges emit 
wavelet that propagate with the sound velocity in all 
the directions and are the only source of the surface 
distortion. In that case the compression stage when the 
nucleus remains spherical lasts until the contact point 
2Rwy — vee © 
be found from the geometry of the system) is higher 


velocity v, = Uy(R— vt) / l (can 
than the sound velocity v, inside the nucleus (Fig.1). 
That first compression stage is characterized by the 
uniform pressure along z-axis and the pressure near the 
contact edge being higher than the pressure in the 
central part?!) After that lateral jetting occurs and the 
pressure gradually decreases while the form of the 
nuclei is changing (Fig. 1(a)). 

In all the calculations a simple approximate 
relation between the sound velocity v, and the particle 
velocity v, is used: 


Us = Ugg + kt, (1) 


S 
where k is a coefficient that is used as an adjusting 


parameter in the model and vso is the ambient sound 


velocity. 
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Fig.1 Initial stage of the drop collision with the rigid wall. (a) 
the picture of the system at 1>tje.; compressed zone (according 
to the impact theory?” is shown in grey; (b) Huygens 


construction shows the shock envelope in the colliding system. 


The compression part of the nuclear matter 
Eos?! 


U(p) — 42 +B 
Po 
2 OE, 


Op 


, 

2 

Po (2) 
Pra g ++ f'vty 

> Lp 5 F 2 p)p 


pg 


together with the continuity equations and momentum 


and energy conservation 


Ci + div(pv) =0 


Ot 

O(pv,) B O(pu,v, + pó;.) Q) 
Ot Oz, 

O(pE) 


Ss + div(v(E + p)) = 0 


are used for phenomena analysis. Parameters A, B and 
y are chosen to give the different values of 
incompressibility? ^l. U(p), Ep, p, P and v are the 
Skyrme type mean field potential, binding energy per 
nucleon, density, pressure and particle velocity, 
addition 


considerations are used to find the relations between 


respectively. In simple geometrical 


the normal and tangential components of sound and 
particle velocities in the disturbed and undisturbed 


parts of the system!) 


The introduced acoustic model is applicable 
during the first stage of the collision until there is no 


3132] This fact allows calculating the 


lateral jetting! 
jetting time from the anomalous behavior of the 
contact edge velocity (Fig.2). One may see that U, 
decreases from infinity to some final value and after 
that starts increasing. Such a picture seems to be 
unphysical and the minimum should indicate the 
boundaries of the model applicability and therefore it 


corresponds to the jetting time. 


vD 


c 


0.4 


0.3 
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Fig.2 Contact edge velocity dependence on the sound velocity 
in the compressed part (normalized to the velocity of light). 
When jetting occurs the shockwave starts to 
travel up along the free surface of the nucleus. This 
allows using the model of a jet formation between the 
rigid wall and the shockwave due to the pressure 
difference and the existing particles flow inside the 
nucleus. In this case assuming that the matter in the jet 
is at normal density po it is possible to start with the 
[56] that 
take into account work against the surface forces L 


hydrodynamic equations for the jet parameters 


with the surface tension coefficient o: 


dL, -8zo RdR, (4) 
- t 
L = 3 RiesP rou fojae. (5) 


0 


From Eqs.(1)-(5) it is possible to calculate the main 
characteristics of the system geometry depending on 
the density at contact point Ppouna: 
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2 4 42 
vt 
R.(t) 1 P bound jet , (6) 
j 32 po 


where Rjet and Vjet are the size of the jet (Figs.3 and 4) 
and the particle velocity in the origin of the jet. 

To check the hydrodynamic approach we have 
chosen “*Nb+Nb system with the impact velocity 
V=0.34c that is the same as in the existing BUU 
calculations!" The parameter k=0.77 is adjusted from 
the comparison of jetting time in nuclear fluid 
dynamics calculations for *’Kr+*°Kr at 400 MeV/A Pf! 


and in our model. 


3 Results and discussion 


All the calculations are provided for the two types of 
the EOS. The results are presented in Tables 1 and 2. 


Table 1 Comparison of the main characteristics of the process 
in the suggested hydrodynamic model and BUU calculations!) 


This work BUU!” 
Jetting time, 7 jet, fm/c 3:5 
Shape at Tje AD / D 0.2 
Pressure release time, fm/c 63 40~80 
Maximum density, p/po ~1.76 ~1.5 


Table 2 Geometry of the system in the cases of soft (K=200) 
and stiff (K=400) EOS 


Soft EOS Stiff EOS 
Shape change, AD / D 0.5 0.4 
Geometry, (D + 27) / D 1.1 1.7 
Time from collision, fm/c 24 21 


First thing to mention is the good 
correspondence of the pressure release time (time 
when the expansion shockwaves reach the center of 
the nuclei and therefore the density of the system 
drops down) obtained in our calculations with that in 
the BUU (Table 1). Our calculations have shown that 
from the hydrodynamic point of view cases of “soft” 
and "stiff" EOS differ not only quantitatively but 
different mechanisms are involved in the system 
evolution. The geometry of the system for the two 
cases at a moment of time when the shockwaves reach 
the top of the nuclei and its evolution are shown in 
Figs.3 and 4. From the comparison one may see 


difference in the qualitative pictures for different EOS. 


Fig.3 System geometry in the case of soft EOS. (a) time when 
the shockwaves reach the free surface of the nuclei; (b) cavity 
formation due to rarefaction. Future “bubble” is shown in black. 


Fig.4 System geometry in the case of stiff EOS. (a) time when 
the shockwaves reach the free surface of the nuclei; (b) rim 
formation due to the surface and viscous forces. “Doughnut” 
structure is formed during the system evolution. 


Soft EOS (Fig.3): 

Such geometry allows the shockwaves produced on 
impact to be reflected from the free boundary and 
focused on the symmetry axis of the system. In that 
case rarefaction with cavity formation along the axis 


occurs (Fig.3(b)). Next steps in such a scenario are the 
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cavity collapse and bubble entrapment inside the 
system. This mechanism shows the possibility of the 
bubble formation in the collision. 

Stiff EOS (Fig.4): 

In that case system looks not like a spherical nuclei 
with small side parts (as it is for soft EOS) but rather 
like an expanding “pancake”. In our opinion this 
picture has some similarity with the problem of 
expanding liquid sheet?" Therefore we suggest that in 
this case the rarefaction along the symmetry axis does 
not play that important role but one rather observes the 
“pancake” becoming thinner and the thicker rim being 
formed due to the surface and viscous forces. 

Both mechanisms give low values of the 
expanding velocities at the final state corresponding to 
the low kinetic energies of the fragments obtained in 
the BUU calculations. 

We would also like to mention the possibility 
of the other geometry such as uniform system for 
different values of the impact energy or intermediate 
of the 
occurrence as well defining the energy range where 


values incompressibility coefficient. Its 
the exotic topologies can be observed require some 
further studies of the possible system splashing or 
recoil. Intuitively, the difference between the observed 
qualitative pictures has the origin in the initial stage 
when the compression pressure and temperature are 
different for different incompressibility coefficients. 
At this stage the only governing parameter is the 
sound velocity in the system. Viscous and surface 
forces come to play later on when the geometry has 
been defined already and govern the fragment 
formation processes. 

Therefore, we suggest that the exotic structures 
formation in the head-on heavy ion collisions can be 
explained from the hydrodynamic point of view. In 
order to explain the breakdown of such structures and 
to predict the fragments mass distribution it is 
necessary to consider the possible Rayleigh-Plateau 
instabilities, capillary waves and the influence of the 
long range Coulomb forces at the final stage of the 
system evolution. 


4 Conclusion 


Suggested approach allows for a simple physical 


picture of the exotic structure formation in the head-on 
heavy ion collisions. The qualitative picture obtained 
within our model is in a good correspondence with the 
one observed within the BUU calculations and in the 
experiment. The straightforward link between the EOS 
and the exotic shapes of different topology is 
explained from the hydrodynamic point of view. 
Combination of the introduced hydrodynamic 
approach together with the microscopic calculations 
nature of the 

and give the 


can reveal the physical 
multifragmentation 
possibility to extract the data on the EOS from the 


head-on heavy ion collisions. 


phenomena 
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